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Abstract 

We formulate in terms of the quantum inverse scattering method the exact solution of 
a spl{2\\) invariant vertex model recently introduced in the literature. The corresponding 
transfer matrix is diagonalized by using the algebraic (nested) Bethe ansatz approach. 
The ground state structure is investigated and we argue that a Pokrovsky-Talapov tran- 
sition is favored for certain value of the 4-dimensional spl{2\l) parameter. 
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1 Introduction 



Recently a novel solution of the Yang-Baxter equation invariant under the spl{2\l) symmetry 
has been obtained in the literature by Bracken et al and Maassarani §|. Its peculiar 

feature is the presence of an additional non-additive parameter whose origin goes back to the 
continuous 4-dimensional irreducible representation of the algebra spl{2\l) 0] The physical 
meaning of this solution is that the associated quantum one-dimensional Hamiltonian can 
be interpreted as an exactly solved model of strongly correlated electrons possessing extra fine- 
tuned hopping terms and electron pair interactions besides the typical correlations appearing 
in the Hubbard model. Most of the Bethe ansatz results concerning this model were obtained 
either by using the coordinate Bethe ansatz method directly on the fermionic Hamiltonian ||^, |^ 
or by taking advantage of the fusion construction of the 4-dimensional spl{2\l) R- matrix. 0. 
In this last case, we note that the eigenvalues of the corresponding transfer matrix and the 
Bethe ansatz equations have appeared only as reasonable conjectures. We also remark that the 
authors |p managed to present the phenomenological analytical Bethe ansatz approach for a 
special point of the continuous parameter in which the i?-matrix can be seen as a braid-monoid 
invariant. However, for this interesting system, a more unified approach such as the quantum 
inverse scattering framework [0, |Tl|, |12| has not yet succeeded to be formulated. We remark 
that the quantum inverse scattering approach makes it possible to present the Bethe ansatz 
results in an elegant way: the Bethe ansatz equations and the eigenvalues of the transfer matrix 
appear as a consequence of systematic algebraic manipulations of the creation and annihilation 
operators acting on a certain pseudovacuum. The main purpose of this paper is to fill this 
gap by presenting a detailed formulation of the quantum inverse scattering method for such 
new i?-matrix invariant by the superalgebra spl{2\l). For sake of simplicity, we are going to 
describe our formulation for the rational limit of the sj9/(2|l) i?-matrix. Some numerical and 
analytical results concerning the ground state structure of this model are also presented. 

This paper is organized as follow. The next section is concerned with the presentation of 
"'^We remark that the i?-niatrix in the context of the fundamental representation of OSP{2\2) algebra has 
been previously discussed in ref. H 
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the Boltzmann weights, the associated Hamiltonian and the basic properties of the spl{2\l) 
rational i?-matrix. In section 3 we formulate the quantum inverse scattering method and in 
section 4 we obtain several important commutation relations. In section 5 we elaborate on 
the construction of the eigenstates and it is shown that the eigenvalues depend on an extra 
auxiliary problem of diagonalization. The corresponding nested Bethe ansatz equations and 
the eigenvalues are then determined in section 6. In section 7 we discuss the ground state 
structure of this model and we present some evidence that a Pokrovsky-Talapov transition 
may occur for a special value of the non-additive parameter of the spl{2\l) i?-matrix. Section 
8 is reserved for our conclusions and remarks concerning the application of our formulation to 
other integrable systems. In appendices A, B and C we have collected some useful relations 
concerning the sp/(2|l) algebra, the two-particle state, and the three-particle state, respectively 



2 The spl[2\l) vertex model and its properties 

We start this section by describing the Boltzmann weights of the sp/(2|l) i?-matrix recently 
proposed in refs. ||T], ^ We shall consider the rational limit of the sp/(2|l) i?-matrix, avoiding 
certain extra mathematical manipulations typical of trigonometric weights. We note, however, 
that the vertex structure is basically the same for both rational and trigonometrical cases and 
therefore many of the ideas described in this paper are quite general. Here we will adopt the 
notations of the appendix of ref. The rational spl{2\l) i?-matrix can be written in terms 
of certain combination of projectors Pi{h) and -P3(&) of the sp/(2|l) algebra as follows: 

where A is the spectral parameter, b characterizes the continuous 4-dimensional representation 

of the sp/(2|l) algebra and I is the identity operator. For completeness we present the explicit 
^We remark that at points h = ±1/2 the corresponding one-dimensional Hamihonian is singular. Here we 
exclude of our analysis such atypical representation. 
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matrix expression of the projectors Pi{b) e Psib) in Appendix A. Using these formulas one can 
see that the rational sp/(2|l) i?-matrix consists of 36 non-null Boltzmann weights. We have 
schematized them in figure 1. We choose to represent each bond of the lattice with the variables 
fi and hi {i = 1, 2) in order to represent the fermionic (/i, and the bosonic (bi, 62) degrees 
of freedom used by Maassarani More precisely, the spl{2\l) i?-matrix can be considered as 
a matrix acting on the tensor product of the two 4-dimensional auxiliary space x and 
can be arranged as a 16 x 16 matrix which in the fbbf grading possesses the following form 
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In the expression (2) the functional dependence of the Boltzmann weights is given by 



m(A) 
/(A) 



1 - 26 - 2A 
1 -26 + 2A' 

1-26 
" 1-26 + 2A' 

2A 

1 -26 + 2A' 



r(A) = 
g(A) 



1 + 26 - 2A 
1 + 26 + 2A' 

1 + 26 
" 1 + 26 + 2A' 

2A 

1 + 26 + 2A' 



n(A): 
p(A) 
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;i-26)(l + 26) 
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-4A(A + 1) 
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^ 2A[(1-2&)(1 + 2&)]V^ ^ l-4&^ + 4A 

"^^^ (l + 26 + 2A)(l-26 + 2A)' (1 + 26 + 2A)(1 - 26 + 2A) ^ ^ 

A general feature of the i?-matrix (2) is that many of their Boltzmann weights are not 

invariant under the charge 1^2 symmetry for same specie of bosonic and fermionic index, 

namely /i ^ /2 and 6i ^ 62- For instance, from figure 1 and expressions (3), these are the 

cases of the following pairs of weights {/(A), r(A)},{m(A), q{X)} and {/(A), (7(A)}. In such case, 

a standard crossing symmetry cannot be implemented, since it is not possible to find common 

crossing factor for the weights Z(A) and r(A), for example. In fact, under the 1^2 symmetry 

the parameter b is reflected to —b. Hence, the remaining invariance is the reflection b —b, 

and the physical properties derived either of -R(A, b) or of -R(A, —6) should in fact be the same. 

We remark, however, that the point 6 = is clearly an exception of the discussion we have 

made above. For this special point, the 1^2 symmetry is present, usual crossing relations 

are then possible to be established and consequently a interpretation of -R(A, 6 = 0) in terms 

of a factorizable S'-matrix turns out to be possible 0]. 

Further properties of the i?-matrix (2) can be seen by using certain relations satisfied by 

the projectors -Pi (6) and ^3(6). We have collected such important relations in Appendix A. 

For example, one can show that the quasi-classical r-matrix originated from the i?-matrix (1) 

encodes the spl{2\l) symmetry in a standard way, namely in terms of the spl{2\l) Casimir 

operator. In fact, making the redefinition ^ and expanding around i] = we find 

v{c{b)-iy 



RiX,b,r]) ~ 



1 + 



(4) 



2A 

where r] is the classical parameter, C{b) is the Casimir operator of spl{2\l) and is the graded 
permutation operator. Details of this calculation can be found in Appendix A. Another 
interesting feature of such i?-matrix is its braid-monoid representation at the special point 
6 = 0. In particular, a Temperley-Lieb operator E (monoid) can be defined in terms of the 
projectors -Pi (6) and -P3(6) as follows 

E = lim46[Pi(6) -P3(6)] (5) 

b—*0 

which together with the permutation operator P^ satisfy the braid- monoid relation M. In this 



sense, it seems very interesting to find the underlying algebraic structure for the points b ^ 0. 
After Baxterization such general structure should produce the spl{2\l) i?-matrix (2). A precise 
answer to this question has eluded us so far. 

3 The quantum inverse scattering approach 

The main purpose of this section is to begin the formulation of the eigenvalue problem of the 
corresponding transfer matrix T(A) of the spl{2\l) vertex system on a square lattice of size 
L X L. The diagonalization problem 

r(A)|$) =A(A)|$) (6) 

can be solved by using an algebraic construction [|TU|, [TT|, based on the Yang-Baxter algebra 
of monodromy matrices T{X) 

R{X - fj.)T{\) ® r(/i) = r(/i) ® T{\)R{\ - ^j) (7) 

where the matrix T(A) acts on the tensor product of an auxiliary space and a quantum space 
® and is given in terms of the product of vertex operators C{X) by 

r{x) = Codx)CoL-i{\)....Coi{\) (8) 

where the index 'o' stands for the 4x4 auxiliary space, and as usual the transfer matrix T(A) is 
obtained as a trace of the monodromy matrix T(A) over such auxiliary space. The elements of 
the vertex operator L'^{\) are related to those of the spl{2\l) R- matrix (2) by a permutation 
on the C'^ x tensor space 

LltW = <(A) (9) 
The corresponding quantum one-dimensional Hamiltonian can be obtained as the logarithmic 
derivative of the transfer matrix T(A) at the regular point A = 0. After some algebraic 
manipulation (see Appendix A), the associated spin chain can be only written in terms of the 
Casimir operator of the spl{2\l) algebra by the following expression 

H = - ^2b+mi-2by ^}^^^ + ^^'^^ " + ^b')C.Mb) - Cl^iib)] (10) 
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where we assumed standard periodic boundary conditions. As a consequence, we can see that 
the reflection symmetry 6 — — 6 is exphcitly exhibited by expression (10) since there exists an 
isomorphism between Cji+i(6) and Ci,i+i(— 6) ( see end of appendix A ) . 

Before going on, we remark that the i?-matrix (2) is a null-parity (Grassmann) braid 
operator, and after some sign definitions^ produces a vertex operator which solves the graded 
Yang-Baxter equation |]13|. In this case one has to use the supersymmetric formalism developed 



in refs. [0 by basically changing standard properties such as trace and tensor product by their 
analogs on the graded spaces. However, a graded formulation does not simplify the original 
problem of diagonalization of the transfer matrix. On the contrary, in terms of practical 
calculations with the corresponding Hamiltonian, one has to be very careful to keep track of 
the fermions signs appearing on the tensor product of the Hilbert space. Here we would like 
to stress that the diagonalization of the corresponding Hamiltonian for small lattice sizes is 
extremely important as a guideline, giving us a correct insight of basic properties such as the 
ground state structure. For that reason we stick with the standard formalism, in which such 
job can be performed in a more direct and safe way. Anyhow, the basic difference between the 
standard and graded formulation will be the presence of extra phase factors on the Bethe ansatz 
equations. However, such phase-factors can be accomplished as general boundary conditions 



similarly as has been done before for the graded 0SP{1\2) model |n . 

After presenting the basic definitions and discussions concerning the diagonalization prob- 
lem (6), we are going to turn our attention to the construction of the eigenstates |$) and 
eigenvalues A(A), respectively. We start our discussion by solving the commutation relations 
which follows as a consequence of the Yang-Baxter algebra. 



^The graded _R-niatrix is related to the standard one (2) by R^fiX) — (— l)P^"^P'''^i?^^(A), where denotes 
the Grasmmann parity oi i — 1, - ■ ■ ,4. 
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4 The fundamental commutation relations 



The proper way to work out the intertwining relations (7) depends, to some extent, on the 
properties of the vertex operator JC{X) when it acts on a given reference state. Let us consider 
as the reference state the usual ferromagnetic pseudovacuum given by 

/1\ 



io)=n®io).> io). 



i=l 








(11) 



By using equations (2) and (9) we find that the vertex operator jC{X) acting on the state 
|0)^ has the important triangular form 

f /(A) * * * \ 
/(A) 
/(A) * 
V p(A)/ 
Now, if we write the monodromy matrix T(A) as a 4 x 4 matrix having the following 
particular form 



|0). 



(12) 



r(A) 



f B{X) Si(A) B2iX) F{X)\ 

Ci(A) Au{X) Ai2(A) Ei(A) 

C2(A) A2l(A) A22(A) ^2(A) 

\CsiX) C4(A) C5(A) D{X)J 
the problem of diagonalization of the transfer matrix becomes 

[i?(A)+EA„a(A) + D(A)]|$) = A(A)|$) 



(13) 



(14) 



Moreover, as a consequence of definition (8), we find that the following diagonal relations 
are also satisfied 

S(A)|0) = [i(A)]^|0), D(A)|0)=p(A)^|0), A„„(A)|0) = /(A)^|0),a = l,2 (15) 
as well as the annihilation properties 

Ci{X) |0) = (i = 1, • • • , 5), Aab{X) |0) = (a ^ 6 = 1, 2) (16) 



In particular for the eigenstate |0) the eigenvalue A(A) is determined to be 



A{x) = m]^ + Y.[fW]' + [pWY 



(17) 



a=l 



In order to construct other eigenvalues one has to find the commutation rules between 
the operators appearing in definition (13). Unlikely to what happens to the 6- vertex model 
[p!0| , |lT| |12] and its multi-state generalizations [|15], |l6l, some of the 'nice' commutation relations 
of the spl{2\l) vertex model are more complicated and require additional work. In many 
cases one needs to combine in a special way certain relations in order to get the appropriate 
commutation rule. Let us illustrate the main idea for the particular case of the commutation 
relations between the operators ^^^(A) and Bc{X). As usual, we substitute the form of the 
monodromy matrix (13) in the intertwining equation (7) and by using the Boltzmann weights 
of the i?-matrix (2) we find the relation 



r1^(A - u)B4^i)AUX) 



/(A - /i) 

a{X-fi) 

—^bc- 



m(A — /i) 



Bb{X)Aac{fi) 



{B{fi)EaiX) + Fifi)CaiX)}, a,b,c= 1,2 



7(A-/^) 

where here and in the following repeated indices denote the sum operation. The elements of 
the vector ^ are 

e = (0 1 -1 0) (19) 

and the matrix f(A) has the structure 

/I 0\ 

s(A) t(A) 

t{X) s(A) 

VO 1/ 



f(A) 



(20) 



Previous experience with the utilization of commutation relations for other systems |]I0|, 0, 
\T% , p!5t |16| suggests us that the term B{u)Ea{X) has a wrong order in the commutation rule 
(18). This can be disentangled by commuting the operators -B(/i) and Ea{X) with the help of 
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the following additional commutation relations 



B{fl)Ea{X) 

cr{X- fi) 



EaiX)Bifi) + !!^j^J^FiX)Caifi) 



^bcBb{l^)Aac{X) 



n{X — n) 



p{X - ^) p{X - fi) 

obtaining final result the expression n 



p{X - n) 
F{fi)CaiX), a = 1,2 



(21) 



Aab{X)B,{fx) 



r':,{X - fx)B,{fx)AUX) 



'^^'-^lE.iX)Bi,) + ^^FiX)CM 



m(A — fi) 
1 

p{X - fi) 



5,(A)A,,(/i)-e 



cr(A — fi) 



bc- 



[p{X-fi) PiX-fi) 
where now the matrix r has the following form 

/I 0\ 

6(A) a(A) 



/(A-/i) 
F(/x)a(A)|, a,6,c=l,2 (22) 



r(A) 



a(A) 6(A) 

Vo i; 



(23) 



with 



a(A) 



A 



-, 6(A) 



1 



(24) 



A+r ' ' A+i 

It is remarkable that such r-matrix is precisely the one appearing on the isotropic 6-vertex (or 
the XXX Heisenberg) model. In some sense the Boltzmann weights of the sp/(2|l) i?-matrix 
(2) conspire together in the commutation rules in order to give us as a hidden symmetry 
the 6-vertex structure. We think that this is the algebraic reason why the bare scattering 
matrix of the equivalent fermionic problem (in the coordinate Bethe ansatz approach) |^ 
was previously determined to have the 6-vertex form. The other commutation rules concerning 
the operators Ba{X) which are important in the diagonalization problem (6) can be obtained 
by basically following the method described above. They are given by 



BiX)BM 



^^-5.(/i)5(A) - ^?^^^5a(A)i?(/i), a = 1, 2 



/(/"-A) 



/(/i-A) 



(25) 



*Here we have used the identity n{x) — p{x) — 1. 
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p(A-/x) P(A-/^) 



Ba(X)B,(^i) = y7v^C'(A - /.)Se(/.)S,(A) - a ^^'^ 



Z(A-/.)'-^- ^""/(A-/.)p(A-/.) 

{/(A - fi)F{\)B{fi) - F{fi)B{X)}, a = 6 = 1, 2 (27) 

The commutation rules (22,25-27) form the basis of the algebraic Bethe ansatz for the 
creation operator Ba{\). It turns out, however, that F{X) also works as a creation operator, 
and, as we shall see in next section, it plays an important role in the eigenstate construction. 
Therefore, the commutation relations with the operator F{X) are also necessary and they are 
given by the expressions 

"'^ - '^'-BMEM + '1' " ^> f ' - F(X)AM, a =1,2 (28) 



'^Wn,) = '^/Mm - '^^i^W^) + ^^^MX)EM (30) 

Other necessary commutation relations are mentioned in the appendices B and C. In the 
next section we shall use all of them in the construction of the eigenstates |<l>). 



5 The construction of the eigenstates and the eigenval- 
ues 

We now have almost the complete machinery to start the construction of the eigenstates of 

the sp/(2|l) transfer matrix. The eigenstates can be obtained by acting the creation operators 
Ba{X) and F{X) over the ferromagnetic pseudovacuum |0). Instead of presenting the general 
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solution, we find more illuminating first to discuss our method of construction for the lowest 
eigenstates. The first excitation over the pseudovacuum |0), i.e. the one-particle state, is given 
in terms of a linear combination between the operators Ba{X) by 

|$i(Ai)) =S,(Ai)F'^|0) (31) 

The diagonalization problem (6) for the one-particle eigenstate (31) is solved by making 
use of the commutation rules (22,25-27) and of properties (15,16), and we find the following 
important relations 

B{X) 1$ i(Ai)) = J^^^^im' l<^>i(Ai)) - 7(^3^ [^(^i)]'^^'^^^" |0) (32) 

m 1$ i(AO) = 1[IzI'] ipW]' i'^i(^i)) + ^[^Z^|] [/(Ai)]^ai^a(A)F^ io) (33) 

E ^.(A) |$i(Ai)) = [1 + a(A - Ai)][/(A)]^ |$i(AO) 

a=i n-^-M) 

"^^^ ~ ^^^-[/(AO]"fia(A)F'^ |0) - 4^^[/(A0]"ai5a(A)F'' |0) (34) 



fiX-X,)''' ^ P(A-Ai) 

The terms proportional to the eigenstate |$i(Ai)) are denominated 'wanted terms' and 
contribute for the eigenvalue A(A, Ai). The remaining ones are the so called 'unwanted terms' 
and they must be canceled out. From expressions (32-34) we can see that this is the case, 
provided that f\ 

where we have used the reflection property 

m(A) m(— A) 

The two-particle state |$2(Ai,A2)) depends both of the operators Ba{X) and F{X). This 
becomes clear if we consider the commutation rule (27), suggesting the following ansatz 

|$2(Ai, A2)) = S,(Ai)5,(A2)F^" |0) + /i(Ai, A2)[/(A2)]^F(Ai)aaF'" |0) (37) 



(36) 



^ Here it is interesting to point out that condition (35) does not impose any further restriction on the 
constants appearing in the linear combination (31). As a consequence, the corresponding eigenvalue is 
double degenerated ( see also section (7) ). 
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The state |$2(Ai, A2)) degenerates several unwanted terms if one tries to solve the corre- 
sponding eigenvalue with the help of the commutation rules (22,25-27) and (28-30) for the 
operators Ba{X) and F{\), respectively. There exists some unwanted terms which can be au- 
tomatically canceled out by an appropriate choice of the function h{Xi, A2). For instance, the 
unwanted terms 

are all of them proportional to (F^^ — F^^) and can be excluded by fixing the following form 
for the function h{Xi, X2) 

h{X^, A2) = h{X^ - A2) = -4y^^ (39) 

P[Ai — A2) 

The other remaining unwanted terms need a further restriction in order to be canceled out. 
For example, by collecting the contributions concerning the term of kind Ba{X)Bh{X2) we find 
that 

""^^ ^ [/(Ai)]Mr(Ai - X2)F-^ - [/(A0]^^^^^F^4 B^{X)B,{X2) |0) 



Considering the identity 

/(A) 1 



/(A) /(-A) 

the term given in equation (40) is canceled by imposing the following restriction 



(40) 

(41) 



l{Xi 



L 



f^^^^^ F'"^ = -aK-X,)F"^\ t^j (42) 

In fact, in Appendix B we show that all unwanted terms of many different kinds can be 
canceled out by using such restrictions for both Ai and A2. We remark that conditions (42) arc 
particular cases of the general Bethe ansatz equations which are going to be discussed in the 
next section. Before going on, it is important to notice that a certain recurrence relation can 
be established between the one and the two-particle states. In order to see this, it is convenient 
to write our results in a more compact way. Let us define the n-particle state as 

|^n(Ai, • • • , A„)) = $n(Ai, • • • , A„).F |0) (43) 
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— * — * 

where • ■ ■ , A„) and F are vectors with 2" components. Here we shall denote the com- 

ponents of vector F by F"'"-''"'"'^ . Considering the bilinear vector B{X) as 

5(A) = (5i(A) 52(A)) (44) 

and taking into account our previous results (31,37) and (39), the corresponding two-particle 
vector $2(Ai, A2) is then written as 

<S2(Ai, A2) = 5(A) ® $i(A2) + [/(A2)]^4T^^^F(Ai)f ® $0 (45) 

— * 

where $0 is the unitary constant. A remarkable property present in equation (45) is the 
symmetry under the exchange Ai ^ A2. More precisely one can show that 

$2(Ai, A2) = $2(A2, >^if-^Y^^ (46) 

l[Ai — A2) 



where we have used the commutation rules (27) and the identity 

rll{X)-rll{X) (47) 



K-x) 

The exchange symmetry between the variables A^ is always a welcome feature in the al- 



gebraic Bethe ansatz analysis (see e.g. refs. |jT2|, [T^)- Such property can be used to cancel 



several kinds of unwanted terms differing under the permutation on variables {Aj}. By using 
this property we can avoid extra cumbersome mathematical analysis, like that appearing on 
the direct proof we gave in Appendix B for the variable A2 entering in the two-particle state. 
This last discussion and the recurrence relation (45) serve as a motivation for us go on and 
to search for the three-particle eigenstate. As before, we start with an ansatz which is able 
to collect together the 'easy' unwanted terms. We find that the three-particle state has the 
following structure 

$(Ai, A2, A3) = S(Ai) ® $2(A2, A3) + [/(A2)]^F(Ai)4^^^^e ® $i(A3)F2(A2, A3) 

p{Xi - A2) 

+ [/(A3)]^F(Al)4^^^® $i(A2)F3(A2, A3) (48) 
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A simple way of determining the unknown function Fj{X2, A3) is by combining the exchange 
symmetry property with the direct cancellation of the simplest unwanted terms. From the 
permutation symmetry we find the relations 

m..h) = '-f^^^ (50) 

J[A2 — A3) 

For completeness, in Appendix C we have presented some details of this analysis as well 
as how certain unwanted terms are easily canceled out. Now, after this detailed exposition of 
our method, the general n-particle eigenstate can be, then, obtained as an induction of the 
expressions (45) and (48). As a final result, we obtain that the ra-particle state satisfies the 
following recurrence relation 



4(Ai, ■ ■ ■ , A„) = 5(Ai) ® $„_i(A2, ■ ■ ■ , A„) + F(Ai)f ® 



i=2 



J p(Ai-A,) 



11 Jn rT'I'n-2lA2,---,Aj_i,Aj+i,---,A„j _[J_ — — — — [bi) 

1 k k+l 

where the index under the r-matrix denotes its action position on the tensor product ® ■ ■ ■ ® ® 

n 

■ ■ ■ <S). Under a consecutive permutation Xj ^ Aj+i, the ra-particle eigenstate satisfies the 
symmetry relation 

^n{Xl, ■ ■ ■ , Xj, Aj+i, ■ ■ ■ , A„) = $ri(Ai, • ■ ■ , Aj+i, Xj, ■ ■ ■ , An)^^77T^-^— : ^T"^ (52) 

L[Aj - Aj+iJ 

Here we would like to make some remarks. It is interesting to notice that the recurrence 
relation (51) resembles in some way the one found by Tarasov in the algebraic solution of 
the Izergin-Korepin model {A^ system |18|). Of course, our case is much more involved due 
to the necessity of vectorial notation and the presence of a hidden symmetry enhanced by the 
vector ^ and the 6-vertex r-matrix. Disregarding the obvious subtleties of our construction, 
one is tempted to interpret (51) as the simplest vectorial generalization of the one found by 
Tarasov for the model. In fact, the 6-vertex model is one of the simplest system possessing 
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a non-diagonal and factorizable r-matrix. In this sense, we believe that our recurrence relation 
(51) is an important guideline for further generalizations concerning the presence of other 
non-diagonal r-matrix. 

In order to close this section, we turn to the restriction condition on the variables {Aj}. 
An important unwanted term comes from the structure Baj^{Xi)Ba2{^2) ■ ■ ■ Ba„{Xn) when one 
of the Xi is exchanged with A, producing for example the Ba^{X)Ba2{\2) ■ ■ ■ Ba,X^n) unwanted 
term. Such term is generated by the action of the operator J2a ^aa(A) and -B(A). By using the 
commutation rules (25) the term coming from -B(A) is 

^^^^[^(Ai)]" ft ■ ■ •^-(^") |0) (53) 

and under action of J2a ^aa(A) we get 

- ^^^^^[/(Ai)]" ft f^x,^_xf ^cr.n^' - A,)rS-(Ai - A3) ■ ■ ■ rtjZ.A^i - A„) 

F"^B,^_,{X)B,,{X,) ■ ■ ■i?e„_,(A„) |0) (54) 

One can write this last term in a more compact form by defining an auxiliary transfer 
matrix associated to the problem of an inhomogeneous 6-vertex system as 

T^'\K {A.})r.t = «(A - Ai)<-(A - X,) ■ --rtri^ - K) (55) 
and equation (54) becomes 

- ^^^^^[/(Ai)]" ft ^(^^^_^.^ tW(A = A„ {A,})-:.tF--- -i?,,(A)i?,,(A2) ■ ■ ■ i?,„(A„) |0) 

(56) 

Here we remark that the same kind of reasoning can be done for any Aj, namely if Ba-{^) 
replaces BaX\)- Now one has to use the property (52) and to perform cyclic permutations 
until one gets the variable Aj on the first place. The basic difference is that now we get a string 
of ordered r- matrices multiplying the components F"'"'""'^ of vector F. This argument is com- 
monly used in many other algebraic constructions [JI2|, and it stresses the importance of the 



exchange symmetry (52). Taking into account such discussion and collecting together equa- 
tions (53) and (56), we find that the unwanted term Ba^{Xi) ■ ■ -BaX^) " " '-San (An) is canceled 
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provided that 



L 



n = T(^(A = A., {A,})^;.iF^"-^^ (57) 



which for = 1,2 reproduce our early conditions (35) and (42), respectively. 

Other crucial unwanted term is that made by replacing one Ba-{^) by £'a^(A), such as 
Eaj^{\)Ba2{^2)---Ba„{Xn)- These terms come out the contribution of J2a^aaW and -D(A). The 
piece which comes from J2a^aa{^) has the form 

p{X - Ai) ^.^2 /(^i - ^i) 

and the contribution of -D(A) is 

- ^|^^[/(Ai)]" ft A,) '"^"'^^^ ~ ^2)rt-(Ai - As) ■ (A, - A„) 

F"^^ikEkiX)B,A>^2) ■ ■ ■i?e„_,(A.)5z,,„., |0) (59) 

By combining these two terms together and by using definition (55) one concludes that 
such unwanted terms are canceled out by the same restriction (57), as it should be. Of course 
we have several others non-trivial unwanted terms. Although we do not have a systematic 
proof that such terms vanish, the checks we performed so far have been rather exhaustive ( 
see e.g. appendices B and C ), and seem to us to establish beyond reasonable doubt that the 
restriction (57) is the unic condition to be imposed on the n-particle eigenstate in order to 
cancel the unwanted terms. 

Finally, the eigenvalue A(A, {Aj}) of the ra-particle eigenstate can be calculated by keeping 
only the terms proportional to the eigenstate |$„(Ai, ■ ■ ■ , A„)). For instance, that proportional 
to the vector Ba^{Xi)Ba2{X2) ■ ■ ■ Ba„{Xn) is determined by the extensive use of the first terms 
of the commutation rules (22,25-27), and we find that 

(60) 



^Some trivial unwanted terms come in pairs and are automatically eliminated such as we have exemplified 
for the cases of the two and three-particle states. 
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where A*^^^(A, {Aj}) is the eigenvalues of the auxihary problem related to the transfer matrix 
of the inhomogeneous 6-vertex model, i.e. 

In conclusion, the results of this section show us that the computation of the eigenstates 
and the eigenvalues of the spl {2\1) model is still dependent of an additional eigenvalue problem 
for the inhomogeneous transfer matrix T^^\X, {Aj}). We shall discuss this matter in the next 
section. 



6 The nested Bethe ansatz equations 



The auxiliary problem (61) can still be solved by an algebraic approach since the Yang-Baxter 
algebra (7) is also valid for inhomogeneous transfer matrices. The corresponding monodromy 
matrix has the form 



T(i)(A, {A.}) = 4^j(A - A„)/:ili(A - A„_i) ■ ■ ■ C^^,\\ - Ai 



'(1)/ 



(62) 



where C^^\\) is the vertex operator of the isotropic 6-vertex model and its matrix elements 
on the space x are given by 

/I 0\ 

a(A) 6(A) 

b{\) a(A) 

Vo 1/ 

In order to go on, we basically have to adapt the well known algebraic results for the 



£«(A) 



(63) 



homogeneous 6-vertex model of Faddeev et al |T0|, O, IT2[ in the case of an irregular lattice. 



We recall that this problem has appeared in many different contexts in the literature [T^, |T6 



but for sake of completeness we present the main results. Taking the monodromy matrix as 

A«(A,{A,}) i?«(A,{A,}) 
C«(A,{A,}) D«(A,{A,}) 



(64) 
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and defining the reference state as 

n 

o«) = n®|o(^))^., o« 

we find the following properties 



i=l 







(65) 



A«(A, {A.}) |0«) = |0«) , D«(A, {A.}) |0«) = n «(A - A.) |o«) , C«(A, {A.}) |o«) = 

(66) 



i=l 



The algebraic Bethe ansatz developed by Faddeev et al |T0|, [11], [T2[ says that all eigenvectors 
of the transfer matrix ^(^^(A, {A^}) + D^^\X, {A^}) can be written as 



$w(/zi,---,/z^)) = ni?^'H/^uA,}) o« 



(67) 



i=l 



for some additional restriction on the variables {/ii}. Moreover, the Yang-Baxter algebra (7) for 
the monodromy matrix (62) and the 6-vertex r-matrix (23) yields the following commutation 
relations 



A«(A,{A,})i?«(/i,{A,}) 



^ 5«(/., {A.})A«(A, {\})-^-^l^,B^'\\ {A.})A«(/i, {A.}) 



a(/i — A) 



a{iJ, — A) 



(68) 



D«(A, {A,})5(^)(^, {A,}) = ^T^S^^H^, {A.})I^(^nA, {A.})-^^[^-4i^^'^(A, {A.})/^«(^, {A.}) 

a(A — /ij a(A — ji) 

(69) 

5«(/x,{A,}),5W(A,{A,})] =0 (70) 

By using the commutation relations (68-70) we can carry y4'^^)(A, {Aj}) -|- D'-^-*(A, {Aj}) 
through all B^^^fXi, {Xj}), and we find that the eigenvalue is 



m -1 n m -I 

aW(a, {A.}, {/.,}) = n -7—^ + n «(A - A.) n -73-— r 

j=i - i=i j=i (^[■^ - f^j) 

if the numbers {/ij} satisfy the following system of equations 



(71) 



Yl a{iJ,j - Xi 



i=l 



11 ^ T, J = l,---,"^ 



(72) 
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Hence, by using the expression for A(^)(A, {Aj}, {fjij}) in equation (60), we finally find that 
the eigenvalues of the spl{2\l) model has the following general form 



A(MA.),W) = [/(A)]^ny^jn^^ 



+ l[a{X-Xi) 



n - 

,=1 «(A - i^j) 



A /(A^ - A) 



f=U(A-AO 



(73) 



_f(AO 



(74) 



and the restriction condition (57) for the eigenstates becomes 

^ rn -1 

where we have used the identity given in equation (41). 

The set of coupled equations (72) and (74) are usually denominated nested Bcthe ansatz 
equations. They can be written in a more symmetric way if we perform the transformation 

6+1/2 



A,; 



A, b - 1/2 ^ 
2 2 ' 2 



and afterwords by performing the rescaling 



^ f^j ^ ~ 



(75) 



(76) 



we get 



A,-(&-l/2)^ 
Xi + {h- l/2)i 

Hk- 1^1- 2i 



Xi- Hi + i 



1=1 fJ'k — fJ'l + ^« 



1=1 ~ f^l 

l^k - Aj - i 



k 



:, i = l,---,n 



, m 



(77) 



j=i fJ'k- Xj + i' 

The eigenvalues A(A, {Aj}, {/i/}), after performing the transformations (75) and (76), can 



1 L 



' " A, + {b- l/2)i -2X ^ fii + (3/2 - b)i - 2A 



be written as 

A(A, {A,}, {/.;}) ^^^^2 _ i,)^ + AJ Xi-{b- l/2)i - 2X fj, f^i-{b + l/2)i - 2X 
^ Xi + {b- l/2)i - 2A ^ll - (6 + 5/2)i - 2A 

fXx^-{b + ?,/2)i-2xlX^Jil-{b+l/2)i-2X\ ' [(l/2-6)i + Aj fi, A^ - (6 - l/2)i - 2A 



^ (l/2-6)i-A j''j^ A, + (6 - l/2)i - 2A ^ 



A(i + A) 



[(l/2 + 6)i + A][(l/2-&)i + A] 



n 



Ai + (6-l/2)i-2A 
\ 'Ai-(6 + 3/2)i-2A 



(78) 
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In particular, the eigenenergies E{L) of the corresponding Hamiltonian (10) can be calcu- 
lated by taking the logarithmic derivative of A(A, {fJ'i}) at the regular point A = 0. More 

we find 



precisely, by performing the operation ^ '^'"^('^'^M'Ia'j}) 



A=0 



V 4(&-l/2) 2L 

In order to conclude this section we would like to make the following remarks. First of 
all, our analytical result (78) can be confronted with the conjecture made by Maassarani 
concerning the structure of the eigenvalues of the spl{2\l) model. Here, one first has to disregard 
the presence of certain phase factors, which is the basic difference between the standard and 
the supersymmetric formulation. We conclude that our analytical result for the eigenvalue 
A(A, {Aj}, {fij}) is in agreement with the conjecture made by Maassarani p|, if one takes the 
rational limit of equation (73) of ref. [Q]. Our second remark is concerned with the extension of 
the Bethe ansatz results obtained in this section for a more general classes of twisted boundary 
conditions. These boundary conditions correspond to the introduction of a seam with different 
Boltzmann weights along the infinite direction on the cylinder. Such weights depend on two 



20 



angles 9i and 9'^ and we represent them by the operator C{X) given by 



//(A) 
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where 001^2 = e • The integrabihty is still preserved and the nested Bethe ansatz equations 
(77) are modified in the following way 

'K-{h- 1/2)1^'' 



\i + {h- l/2)i 

m 

n 



2i 



l^l \ — f^l — I 



1, 



, n 



-i(.01'02) 



n 



k 



, m 



1) 



.^y IJ'k - fJ'i + 2i Atfc 
and the eigenenergy equation (78) remains unchanged. One interesting case of twisted bound- 
ary conditions is when 61 = —62 = d/2. In this case the spl{2\l) algebra is preserved and is 
accomplished in terms of the spl{2\l) generators as follows. The odd generators transform as 



and the even generators behave as 



12) 



S. 



± 



Ss^Ss, B^B (83) 
^In the fermionic version Q of the corresponding Hamiltonian (10) this means that the fermions with spins 
up C|(i) and down C|(«) satisfy the boundary conditions C|_j^(l) = e**i'^C|,|(L + 1), c| ^(1)) — e^*^^'^c| + 1) 
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This special boundary condition is quite useful, because by varying the angle until 'd = 27r| 
(in the sector where {L — n) is even) we can study the spectrum of the system in the super- 
symmetric formulation. In this sense, both the standard and the supersymmetric formulation 
may be related by performing in determined sectors of theory a twisted boundary condition. 
This idea has been used with success to study the spectrum of the 0SP{1\2) spin chain . 

The final remark consists of a discussion concerning the reflection symmetry 6 — — 6 in 
terms of our algebraic formulation. As has already been noticed by Maassarani 0, the R- 
matrix structure (2) allowed an additional reference state such as 

/0\ 



i=l 






The triangular form of the vertex operator C{X) is now given by 

/p{\) \ 



^4) 



5) 



* ^(A) 

* ^(A) 
\ * * * r(A) J 

which can be related to the earlier triangular property (12) by b —>■ —b and /i <-> f2- Such 
relation becomes even more rigorous if one looks for the commutation rules comming from the 
Yang-Baxter algebra (7). In this case, we find that it is more convenient to start with the 
following matrix form of the monodromy T(A) 

fD{X) C^iX) C,{X) C'3(A)\ 
E,i\) i22(A) i2i(A) C^iX) 
E^iX) AuiX) in(A) C'i(A) 
VF(A) B2{X) fii(A) B{X) J 

^Notice that this corresponds to a periodic (anti-periodic) boundary condition for the bosonic (fermionic) 
degrees of freedom. 



r(A) 



(86) 
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By using the Yang-Baxter algebra with this structure, the commutation rules of the diagonal 
terms with operator Ba{X), for example, are given by 

g{\-fi) 9{X-n) 
""^^ ~ ^b{X)C(nU, + ^7^^Um)Aca(n) (89) 



(87) 



p{X - ii) p{X - /i) 

Such commutation relations are the same of that given in equations (22,25,26) by perform- 
ing the reflection symmetry 6 — > — 6 and a — > —a. We have checked that similar situation 
appears for the other commutations rules. Since a —a can be considered as a canonical 
transformation, we conclude that the eigenvalues parametrized by equations (77) and (78) ei- 
ther with +b or with —b, should produce the same spectrum for the sp/(2|l) model. From the 
Bethe ansatz point of view, one would think that this sounds a bit strange, since the topology 
of the nested Bethe ansatz equations are very different for ±6. The correct interpretation is as 
follows. They in fact can produce the same eigenvalues, but of course with different structure 
of zeros. Actually, we have checked this fact by solving numerically the nested Bethe ansatz 
equation for some values of ±6 with a finite lattice size L. 



7 The ground state structure 

The purpose of this section is to study both numerically and analytically the behaviour of the 
ground state of the sp/(2|l) Hamiltonian (10). Due to the reflection symmetry b — >■ —b, it is 
enough to investigate the regimes < 6 < 1/2 and 1/2 < 6 < oo. We begin our analysis by 
investigating the behaviour around the reference states with ferromagnetic feature. In order 
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to keep our previous notation (section 2) we define the four basic states acting on the i-th site 
of the one-dimensional lattice by 



n\ 




/0\ 








/o\ 







1 
































1 







i 




i 




i 





(90) 



The low-lying excitations over the ferromagnetic state 

^'^) = l/i(l)---/i(i^)) (91) 

can be constructed if one takes the linear combination of the states made by replacing one 
state fi{k) by 61,2 (/s), namely 

L 



0{:2) = E«(^) 1/1(1) •••^2(A;)---/i(L)) 



k=l 



While the eigenvalue of the ferromagnetic state is trivially determined to be 

2L 



1/2 -h 



that of the low-lying excitation 0{ 2 ) is related to the solution of the equation 



Ei]^a{k) 



^ -{L-l)a{k)+ ^ 



-[a{k + I) + a{k - I)] 



1/2-6' ' ' ' 1/2-6' 

This relation can be solved by taking a{k) ~ Ae^^^ and the dispersion relation is 

2L 4 

El\(p) = + — cos^(p/2) 

i,2\F) 1/2-6 1/2-6 ' 

where under periodic boundary condition the momenta p assumes the values 

27rn 



-; n = 1, • • • , L — 1 



Analogously, by repeating the same reasoning to the other ferromagnetic state 



(92) 



(93) 



(94) 



(95) 



(96) 



4' 



1/2(1) •••/2(i^)) 
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(97) 



we find the following results 

In the regime l/2<6<oowe have El^ip) < and El^ih) > 0. Moreover, the sign in 
the dispersion relation of E(^2{p){E(12{p)) is positive (negative). This suggests that the ground 
state is given by Eq^ (6) and the excitations over it grows until reaching the upper- bound (6). 
Remarkably enough, the statement concerning the ground state can be put in a more rigorous 



ground. In fact, following Bader and Schilling ||2^ the ground state satisfies the relation 



Eo > LE^ (100) 

where E^ is the lowest eigenvalue of the two-body Hamiltonian In the regime < 6 < 

1/2 we find that 

Eo = (101) 

On the other hand, the expectation value of the sp/(2|l) Hamiltonian in the ferromagnetic 
state 1/2(1) ■ ■ ■ f2{L)) is —j^^, and by the variational principle we have 

Eo < -, r 102 

°- 1/2 + 6 ^ ^ 

By comparing equations (100) and (102) we conclude that in the regime 1/2 < 6 < 00 the 
ground state energy is El^{h). Actually, by analysing the dispersions relation E(\{p) around 
p ~ TT we find that it is three-fold degenerated. The rest of the statement, i.e. that the spectrum 
satisfies —jj^^ ^ E/L < — y/I^' been confirmed by a numerical diagonalization of the 
Hamiltonian (10) up to L < 12. 

In the regime < 6 < 1/2 rigorous results become more involved since the Hamiltonian 
(10) is not hermitian f\. For instance, such 'nice' argument of Bader and Schilling |^ needs 



further elaboration. In this case, the ferromagnetic state (po^ ) has lower energy than the state 



Nevertheless, we have numericaUy checked up to i < 12 that the ground state is real. 
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'q^^ and the sign on the dispersion relation of E(^2{p) is now positive. This may suggest 
that there exist states with much lower energy. Indeed, for 6 = we have found previously 
that the ground state has an antiferromagnetic structure and in the thermodynamic limit the 
energy per site e^o was determined to be Cqo = —8ln{2) ~ —5.5452 < i?g^'-^^(6 = 0) = —4.. 
Furthermore, the numerical diagonalization of Hamiltonian (10) for 6 7^ (0 < 6 < 1/2) shows 
that the ground state jumps into the many possible U{1) sectors of the theory if one varies 
the lattice size L. Let us illustrate this from the Bethe ansatz point of view. As an example, 
let us assume that for sufficient small b some states have the zeros structure similar to that 
previously found for the ground state at b = [0, namely 

A, =e,+^ + 0(e-'^^), fij=^ (103) 

By substituting such structure of zeros in the nested Bethe ansatz equations, (see ref. 0) 
we find that the density of variables p{Oy thermodynamic limit, satisfies the following 

integral equation 

/+00 
^'(e - u)p{u)du (104) 
-00 

where the prime symbol stands for the derivative and the functions and ^p{^) are given 
by 

^-(0 = 2[arctan(^y|-^) - arctan(^^^)], <^(0 = 2 arctan(e/2) (105) 
This integral equation is solved by Fourier techniques and we find that 

^^^^ cosh(7rO + cos[(l/2 - b)7i] ^ ' 

Since the sp/(2|l) Hamiltonian can be interpreted as spin-3/2 chain, the magnetization per 
site is then given by 

where the integers n and m (see Bethe ansatz equation) are the number of variables Xj and 
fii, respectively. From equation (106) we find that -^ = 2(1/2 — 6) and ^ = (1/2 — 6), and as 
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consequence of equation (107) we have|^ 

^ = 36 (108) 

It is also possible to verify that the energy per particle associated with the structure of 
zeros (103) is lower than £'o^(6). In summary, we believe that our results lead to the following 
picture. Strictly at 6 = the system has zero magnetization and presents an antiferromag- 
netic behaviour. As soon as we turn 6 > 0, the ground state gets a finite magnetization and 
is partially ferromagnetic ordered. We believe that this picture remains in the whole regime 
< 6 < 1/2. After the singular point 6 = 1/2(6 > 1/2), the system is then fully ordered in the 
ferromagnetic state 0o^^. In this sense the parameter b plays the role of the incommensura- 
bility such as a chemical potential or a magnetic field indicating that at h = 1/2 the system 
presents a phase-transition of Pokrovsky-Talapov type In fact, the typical quadratic form 
of the dispersion relation appears in E{\{h,p) for the low-lying excitation around p n. It 
seems interesting to understand this picture in terms of the particles and anti-particles of a 
factorizable ^-matrix. We suspect that the point 6=1/2 behaves as the threshold for the 
mass of the solitons and the anti-solitons present in the theory. 



8 Conclusions 

We have shown that the one parameter family of an integrable spl{2\l) vertex model is exactly 

solved by the algebraic Bethe ansatz. The eigenstates have been formulated in terms of the 

creation operators through the recurrence relation (51). The eigenvalues of the corresponding 

transfer matrix is computed by solving an auxiliary problem related to an inhomogeneous 6- 

vertex model. We have discussed how the reflection symmetry 6 — — 6 can be encoded in 

terms of the commutation rules. The ground state picture has been discussed and we have 

presented arguments that the system has a commensurable/incommensurable phase transition 

of Pokrovski-Talapov type. 

It is remarkable that although we have considered b << 1/2, the limit b ^ 1/2 in equation (108) reproduces 
the ferromagnetic structure found for b > 1/2 . 
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We believe that the formulation described in this paper is by no means only particular to 
the isotropic sp/(2|l) vertex model. As we have already commented, the whole construction can 
be generalized for the anisotropic model (trigonometric case) almost directly. The study of the 
phase-diagram of the anisotropic model and in particular the Pokrovski-Talapov transition, 
seems to us to be a very interesting problem. Besides such direct generalization, we have 
reasons to think that our formulation is the cornerstone to solve, by the algebraic Bethe ansatz 
approach, certain integrable models related with the symmetry C„. This should be the case of 



the isotropic Sp{2n) and OSP{2\2n — 2) vertex models p2[. Other model that is still waiting 
for an algebraic solution is the Hubbard model. Since this system also possesses a hidden 
6- vertex symmetry, we strongly believe that our formalism can be adapted in order to give 



the algebraic Bethe ansatz solution of the Hubbard model [g^. We hope to report on these 
problems in a future publications. 
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Appendix A : The spl{2\l) algebra and the R-matrix properties 

The algebra spl{2\l) [|| consists of four even generators {S±, S3, B} and four odd generators 
{V±, V±} ■ Following refs. [§, |] they satisfy the commutation rules 

[Ss,Si] = ±Si, [S+,S^] = 2Ss, {r±,F±} = ±l/2^± (A.l) 

{V+,V^} = -1/2P+, {F+,V_} = -1/2P_, [P±,F±] = ±F±, [P^,V±]=±V± (A.2) 

[P±,V±] = [Pt,V±] = 0; {Vi,Vj} = {Vi,Vj} = 0,z,j = ± (A.3) 

where the symbols [, ] and {, } denote the comutator and the anti-commutator, respectively 
. We also have the identity Pj- = S3 ± B . The 4-dimensional representation 0, ^ of these 
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generators in the bffb grading possesses the following matrix representations 



/ 



\ 



e 



a 

^ 

7 ^ 

/? 



\^ y 
/ 



^ ^ 



-a 

e y 



V 



1/2 -b 
1/2-6 



s+ = 

where 4a7 = 
generators as 



1 ^ 




1 + 26 and 4/3e 



^ ^ 

-/? 



\ 7 y 



5_ 
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1/2 + 6 




/ 












-1/2 




/ 





-1/2 + 6 








\ 



1/2 








































1/2 + 6 


\ 







-1/2 j 








-1/2-6 
(A.4) 






-1/2 + 6 ) 
(A.5) 



(A.6) 



1 — 26 . The Casimir operator is written in terms of these 



a,,+i(6) = 2{S+ ®S^ + S^® S+} + 2{P+ ® P_ + P+ ® P_} + 
4{V1 ®V+ + V.®V+-V+®V^-V^®V^} + Ab'^I 



(A.7) 



In equation (A.7) the symbol ® stands for the supertensor product. More precisely the 

s 

elements oi A® B are 



{A ® BY' 



ab 



■^y{i)p{j)+p{a)p{b)+p{i)p{B) j^^^j^^^ 



(A.8) 



where p{f) is the Grassmann parity of the object / . 

As it has been observed by Maassarani the projectors Pi,i = 1,2,3 of the sp/(2|l) 
algebra play a fundamental role in the construction of the P-matrix . In the isotropic limit, 
following ref. H, we find that the 16 x 16 matrices expressions for the projectors Pj in the 
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fbbf grading are 



Piib) = 



/I 
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P3{b) = 



/O 
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where P2{b) = / - (Pi(6) + Psib)) § and a{b) = (26 + l)/8b, (3{b) = (26 - l)/86 and 7(6) = 
Vl — 46^/86 . Besides the usual projectors identities, Pi{b)Pj{b) = 6ijPj{b), the operators 
Pi (6) and ^3(6) satisfy the following useful relations 

P^ = I-2[P,{b)+P,{b)], P^Pr{b) = -P^{b), P^Ps{b) = -Ps{b) (A.ll) 

where (P^Y^ = (— l)^*-*''^*-"''*5i,A;^j,« defines the graded permutation operator. By using definitions 
(A.7,A.8) and expressions (A. 9, A. 10), the Casimir invariant can be connected to the projectors 
Pi (6) and P3(6) as 

a,,+i(6) = / - 2[Pi(6) + P3(6)] + 46[Pi(6) - P3(6)] (A.12) 

Such relations are important in the study of the classical limit of the spl{2\l) P-matrix . 
In fact, by introducing the quasi-classical parameter 77 in equation (1) as A ^ ^ we have 

R(X,b,7i) = I - ^ -Pi(6)-- ^ -P3(6) (A.13) 

^ ' ' (l-26)r/ + 2A ^ ' {l + 2b)r] + 2X ' ^ ' 
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By expanding this last expression around 77 = we find 

R{X, 6, = 0) = / - 2[Pi(6) + Psib)] = (A.14) 

and 

^^^^'"^^ 1^=0 = \ {Piib) + Psib) - 2b[P^{b) - Ps{b)]} (A.15) 

As a consequence, we find that equation (4) of section 2 foUows from the expressions 
(A. 14, A. 15) and the identity (A. 12). Analogously, the corresponding Hamiltonian can be also 
written in terms of the Casimir operator . Considering that the two-body Hamiltonian 
is determined as the derivative of the i?-matrix at the regular point A = , we obtain 

Now if we consider the identities (A.11,A.12) the square of the Casimir operator is 

CfMb) = I + 8b[Ps{b) - Pi(6)] + 16b'[Ps{b) + Pi(6)] (A.17) 
and now by solving equations (A. 12) and (A.17) for the operators Pi{b) and ^3(6) we find 

= -T7 7T7 rT[(4& + 1)/ - 46Ci i+i(&) - Cf,,Jb)] (A.18) 

1 + 26 166(1 - 26) (1 + 26)^^ ' *,*+iv ;j v ; 

^ [Ub _ 1)/ _ 4bCii+i(b) + Cf., ,(6)1 (A.19) 

1 - 26 166(1 - 26)(1 + 26)^^ ^ '^^^ Ji v j 

These last relations are then used in equation (A. 16) in order to reproduce the Hamiltonian 
expression (10) presented in section 3. 

Finally, we remark that the sp/(2|l) algebra is invariant under the following isomorphic 
transformation 

S±^S±, Ss^Sa, P±^P^, V±^V±, V±^V± (A.20) 

which in terms of the parameter 6 means the refiection symmetry 6 — * —6. This can be easily 
seen from the matrix representations (A.4,A.5,A.6) of the sp/(2|l) generators, provided we also 
perform the canonical transformation /i <-> /2. 
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Appendix B : The two-particle state 

The main purpose of this Appendix is to give extra details that the two-particle state (45) 
we have constructed is in fact an eigenstate under the Bethe ansatz restriction (42). In order 
to collect the unwanted terms we have to use the annihilation property (16) and the following 
additional commutation rules 

Ca{X)B,{i^) = B,{u)Ca{X) - '^^^^[B{X)Aab{i^) - B{i^)Aah{X)i a = 1,2 (B.l) 

C„+3(A)56(/^) = '^(^^BMCb+sW + ^-^^^f^^^)CaU^) 

-ea.^|^^{i^(/x)C3(A) + B{^^)D{X)} - ^|^^^a(A)a+3(/x) 

+ '^^^^iimAia{X)AM, a = 6= 1,2 (B.2) 

' '^^^"^^/x)A„,(A), a = 6 = l,2 (B.3) 



/(A-/X)- 

By combining such commutation rules together with those mentioned in section 4 and 
property (15) we find that the non-trivial unwanted terms have the following structures 

S„(A)S,(A,); E,{X)B,{Xiy, F{X){F'^ - F'') (B.4) 

For the first two cases in equation (B.4) it is enough to fix Aj = Ai, since the other possibility 
(Aj = A2) has already been discussed in the main text (see section 5). We now summarize the 
functional forms which multiply the terms (B.4), and discuss how they are canceled out. 

1. The Ba{X)Ba{Xi) term, for a given index a, appears as 



r m(Ai-A)m(A2-Ai) _ m(A2 - A) C(Ai - A) 
1 /(Ai - A) /(A2 - Ai) /(A2 - A) /(Ai - A) 



F'^"[/(A2)]^ + 



r m(A-Ai)m(Ai-A2) _ m(A - A^) rr(A - Ai)C(A - A^) ] ^ 
1 /(A - AO /(Ai - A2) /(A-A2) /(A-AO/(Ai-A) / 
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and by using the identities (36) and (41), we find 

f m(A,-A) + - ^) '"(^^ - {p(A,)]^ + [/(A.)]^) f- (B.6) 



I /(A2 - A)/(Ai - A) /(Ai - A) /(A2 - A,) 
which is null by taking a = 6 in the Bethe ansatz equation (42). 

2. The Ba{X)Bij{Xi) term for a 7^ 6 is more involved. Here the cases Bi{\)B2{X2) and 
B2{X)Bi{X2) are equivalent. The functional structure which comes from the creation operator 
F(Ai) of the two-particle eigenstate is 

(7(Ai - A)ct(Ai - Aa).,., .,i.„i2 ^21 



p(Ai - A) p(Ai - A2) 

and that coming from Ba{Xi)Bb{X2) are 



[1{X2T{F''-F'') (B.7) 



1121 



r m(Ai-A)m(A2-Ai) _ m(A2 - A) [rl|{X^ - X)F'' + rfi(Ai - A)F^ j . ^ 
1 /(Ai - A) /(A2 - Ai) /(A2 - A) /(Ai - A) I U 2JJ + 



m(A — Ai) m(Ai — A2) m{X — A2) 1 



/(A - Ai) /(Ai - A2) fix - A2) /(A - Ai)/(Ai - A) 



^''[/(A2)]^ 

(B.8) 



The simplest way to see that this term is null is to use the Bethe ansatz equations (42) for 
Aj ~ X2 in order to have only terms proportional to F^^[/(A2)]^ and F^^[/(A2)]^. In particular 
we have 

[l{X2r{F'' - F-') = -[/(A2)]^ {r,^2^(A2 - A,) - r'^{X2 - X,)} (F^^ - F^^) (B.9) 

By making such manipulations and collecting the term proportional to ^^^[/(Aa)]^ of ex- 
pressions (B.7) and (B.8) we have 

[ m(Ai - A) m(A2 - Ai) i2/x ^(Ai - A) o-(Ai - A2) . 12/ x xn 21, 



fix, - A) /(A2 - AO '-(^^ - - KA,-A)p(A,-A2)['-^^^ " - '-^^^ " ^'^^ 
^ /(A2-A)/(Ai-A) t^^"^^^ - A)r2^?(A2 - A,) + r-(A, - A)r-(A2 - A,)]} ^-[/(Aa)]^ 



(B.IO) 



which is in fact null if one uses the identity (47) and the factorization relation 

a(Ai-A)ff(A2-Ai) ^^ , _ m(Ai - A) //;(A2 - AQ ^^ 

- piX, - A) p(A2 - AO - + 1) - /(A, - A) /(A2 - AO ~ 
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m(A2 - A) 1 (A2 - A) 

/(A2-A) /(Ai-A) (Ai-A + 1) 

(B.ll) 



Analogously, the terms proportional to F^^[/(A2)]^ are 

m(Ai - A) m(A2 - Ai) 2i/> > x , m(A - Ai) m(Ai - A2) m(A2 - A) 1 

-r 21 {A2 — M) + -TTT -TTT + 



/(Ai - A) /(A2 - A,) ' /(A - Ai) /(Ai - A2) /(A2 - A) /(Ai - A) 

IrlUX. - A)4'(A. - AO +---(A. - ^ - A.)] - f^l'^, _ 

a(Ai - A) a(Ai - A2) 



-[rlliX^ - A) - r,^?(A2 - A^)] F''[f{X,)]^B.12) 



p(Xi - A) p(Ai - A2) 
Now, by noticing that the first two terms in (B.12) can be simplified as 

m(Ai-A)m(A2-Ai) ow. . 
/(A, -A) /(A2-A0'^^^^^"^^^ ^^-''^ 

one can proceed as we did before. If we use the simphfication (B.ll) and the last relation 
(B.13) we finally find that the term proportional to F^^[/(A2)]^ is also null. 
3. The Ea{X)Bb{X2) term for a 7^ 6 appears as 

[ ct(A - A2) /(A - Ai) ct(A - Ai) m(Ai - A2) \ paa^j^^^^-^i^i ^ 



piX - A2) piX - Ai) p(A - Ai) /(Ai - A2) 
/ a{X - Ai) m(A2 - Ai) a{X - A2) rll{X - X^] 



Ip(A - Ai) /(A2 - Ai) p{X - A2) ^(A - Ai) 
and by using the identity 



i^"1/(A2)]'^ (B.14) 



(B.15) 



([/(A2)]^ + [/(A2)]^) (B.16) 



we are able to simplify equation (B.14) as 

[ a(A-A2)/(A-Ai) a(A-Ai)r»(A2-Ai) ' 
\p{X - A2) p{X - Ai) ^ p{X - Ai) /(A2 - Ai) 

which is automatically null by taking a = 6 in the Bethe ansatz equations (42). 

4. The £^a(A)-Ba(A2) term has a contribution from both creation operators: Bi{Xi)Bj{X2) 
and F{Xi). The contribution coming from F{Xi) is 

g(A - Ai) a(Ai - A2).,,- .,^.^^12 tj^21 



g{X - Ai) p{Xi - A2) 
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[/(A2)]^(F^' - F'^) (B.17) 



and those comming from Bi{Xi)Bj{X2) are 

Ip(A-Ai) /(A2-A1) p(A-A2) ^(A-Ai) f^'^''''^ ^ 

\ a{X - A2) /(A - Ai) (7(A - Ai) m(Ai - A2) | p2i^j^^^^-^L 



p{X - A2) p(A - Ai) p(A - Ai) /(Ai - A2) 



(B.18) 



In order to show that all these terms together are in fact null, one can follow the same steps 
of the procedure which we have used for the term Bi[Xi)Bj[X2) (i 7^ j). However, the crucial 
factorization relation here is a bit different, namely 

a{X - Ai) m(A2 - Ai) 21/, _\ \ - - -^1) ^"(-^2 - Ai) 
p( A - AO /( A2 - Ai) - g{X - Ai) p( A2 - Ai) 

, (t{X - A2) 1 (A - A2 + 1) 



p(A - A2) ^(A - Ai) (A2 - Ai + 1)(A - Ai + 1) 



(B.19) 



5. To collect all the non-trivial unwanted terms proportional to F(A)(F^^ — F^^) is a 
very cumbersome job. The main reason is that all the diagonal operators Ylii=\^aa{^)-iB[X) 
and D{X) give non-trivial contributions which are proportional to the many combinations of 
[^(Ai)]n(- • ^[/(AO]^ + (• • •)[^(A2)]^} and [f{X,)]^{{- ■ •)[/(A2)]^ + (• • •)[^(A2)]n. For instance, 
the terms proportional to [/(Ai)]^ are 

[ o-(A - Ai) m(A - A2) 1 (t(A - Ai) m(A - Ai) m(Ai - A2) \ ..^ 

\ p(A-Ai)/(A-A2)/(A-Ai)/(Ai-A) + p(A-Ai)/(A-Ai) /(Ai-A2)r^^ ^ 



^(Ai — A2) n{Xi — A) (y{Xi — A) m(A2 — A) 



.p(Ai - A2) p(Ai - A) p(Ai - A) /(A2 - A)/(Ai - A) 



[/(A2)]^ 

(B.20) 



while those proportional to [/(Ai)]^ are 

/ a(A- A2)m(A- Ai) a(Ai - A2) n(A - Ai) | 
I P(A-A2) p(A-Ai) ^p(Ai-A2)p(A-Ai) j ^ 
[ m(A-Ai) 1 a(A-A2) 21.^ . ^l 
1/(A-A0 ,(A-A0p(A-A2)[' + ^-^^ - 
a{X - Ai) m(A - AQ m(Ai - A2) _ a(Ai - A2) m(A - Ai) g(A - AQ \ ^ 
p(A - Ai) p(A - Ai) /(Ai - A2) p(Ai - A2) ^(A - Ai) /(A - Ai) / ^ ^ ^ " ^ 
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In order to cancel these expressions one has to use the Bethe ansatz identity (42) for both 
Ai and A2. For instance, we first use such relation for A2 in equations (B.20,B.21), and as 
a result we only have terms proportional to [/(A2)]^{(- ■ + (' ' ^[^(^i)]^}- Now, we 

perform the same operation for Ai in order to eliminate the [/(Ai)]'^ terms. The final expression 
is very complicated, but we have checked that it is null with the helping of the Mathematical'^ 
software. 

This finishes our analysis concerning the unwanted terms, and the results of this Appendix 
together with those of section 5 show that the Bethe ansatz equations (42) are sufficient 
conditions to cancel all of them out. 

For sake of completeness let us also discuss the wanted terms. They are responsible for the 
eigenvalues of the transfer matrix . The wanted terms are constituted by two kinds of creation 
operators and from the two-particle eigenstate expression (45) we have 

Ba{Xi)B,{X2)F'^, /(A2)^4y^^F(Ai) (B.22) 

p{Xi - A2) 

The contribution proportional to the first term of equation (B.22) can be obtained directly 
from the commutation relations (22,25-27), and it is easy to get 

A(A, (A.,) ^ HA)^ n ^ ^ /(A)^ n + p(A)^ n (B.23) 

where A(-^'(A, {A,}) is the eigenvalue of the inhomogeneous 6-vertex system with two sites (see 
equation (61) of section 5 ) . Of course the contribution coming from the second term has to 
be precisely the same as that we have got in equation (B.23) . Such calculation is a bit more 
elaborated , since the action of the diagonal operators J^a'^aaiX) , -B(A) and -D(A) on the term 
Ba{Xi)Bf){X2)F''°- can produce many terms of type F(Ai)(F^^ — F^^) . The analysis is as follows 
. The terms coming from D{X) are 

^g(A-Ai)a(A- A2) r(A - Ai) a(Ai - A2)' 



piXnxnF'^ F'') ^^^^ _ _ ^^^^ _ 

and if we use the factorization identity 

_ g(A - Ai)o-(A - A2)p(Ai - A2; 

- ^2) i=i 



(B.24) 



+ r (A - Ai)p(A - A2) = n 5(A - A,) (B.25) 
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we find that the operator D{X) contributes to the eigenvalue with 

Those coming form B{X) are 

a{Xi - X)l{Xi - X)m{X2 - A) l{Xi - A)a(Ai - A2) 



(B.26) 



p(Ai - A)/(Ai - A)/(A2 - A) p(Ai - A)p(Ai - A2) 



(B.27) 



and by using the relation 

cr(Ai - A)m(A2 - A)p(Ai - A2) 



a(Ai - A2) 
we find that the contribution of i3(A) is 



+ /(Ai - A)/(A2 - A) = p(Ai - A)Z(A2 - A) 



(B.28) 



(B.29) 



=-^1 f(\ - A) 

Finally, the operator J2a^aa{X) generates the terms which carry the hidden 6- vertex struc- 
ture . They have been collected as 



/(A)-/(A2)^(F- - F-) - (¥^^[1 ''^^ " 



]+ 



p(Ai-A2)i/(A-Ai)^ ^^(A-Ai) 
a{X - A2)g(A - AiMAi - X2)[rll{X - AQ + rf^(A - AQ] 
p(A - A2)/(A - Ai)(/(A - Ai)a(Ai - A2) 

cr(A - Ai)m(A - A2)p(Ai - A2) 



(B.30) 



/(A - Ai)p(A - Ai)/(A - A2)a(Ai - A2) j 

Remarkably enough the terms on the bracket of equation (B.30) can be factorized only in 
terms of the 6- vertex r-matrix as 

n Jlx^) - - + - - - - Ai)r,^2'(A - A2)} 

(B.31) 

However, by using the auxiliary eigenvalue equation (61) for two sites, one is able to estab- 
lish the following identity 

A«(A, {X,}){F'' - F'') = {rlliX - X,)rll{X - A2) + rll{X - X,)rll{X - X,) 

-rIliX - Xi)rll{X - A2)} (F^2 _ ^21^ ^3 32) 
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and as a consequence of the expressions (B.31) and (B.32) we have that the eigenvalue contri- 
bution of Ea^aa(A) iS 

Hence, this completes the proof that the two-particle state (45) is the eigenstate of the 
spl{2\l) model . 

Appendix C : The three particle state 

This appendix is mainly concerned with the symmetrization property (52) of the three 
particle state . We begin our discussion with the A2 ^ A3 permutation . In this case, the three 
particle vector (see equation (48) of section 5 ) becomes 

*3(Ai, A2, A3) = ^(Ai) ® B{X2) ^(A3) + [/(A3)]^4t^^^^(^i) ® ^(^3)f 

P[M - A3) 



Cj(Ai - A3) DM M X ^ , r;/x mL<^(Ai - A2 



2 



+ / , i H>^iK ® B{X2)F,{X,, A2) + [l{X2)r / . -^(AOe ® i?(A3)F3(A3, A 

P{Ai — A3J p[Ai — A2) 

(C.l) 

In order to relate this vector with the symmetric one $3(Ai, A2, A3), we turn B{X2) over 
-B(A3) with the help of the commutation rule (27) and as a result we get 

S(Ai) ® B{X2) ® BiXs) = {^(Ai) (8) B{Xs) ® B{X2)+ 

— * 

The last term of this identity cancels out the second term in the vector $3(Ai, A2, A3) and 
the symmetrization rule 

$";(Ai,A2,A3) =$3(Ai,A3,A2).^^^f^^^— ^ (C.3) 

l[A2 — As) 

is valid provided the functions F2(A2, A3) and F3(A2, A3) satisfy the following equations 

^3(A3, A2) = F2(A2, A3)Z(A2 - A3)r2-/(A2 - A3) (C.4) 

39 



hiXs, A2) = F3(A2, X3)l{X2 - X3)r2i{X2 - A3) (C.5) 
but in fact they are equivalent since we have the following inversion properties 

l{x)l{-x) = l, T23{x)t2z{-x) = (C.6) 

Similar reasoning can be implemented for the permutation Ai <-> A2 . In this case however, 
besides turning -B(Ai) over -B(A2), we also have to turn -B(Ai) over F(A2)C- By using the 
following commutations rule 

This leads us with the term 

l{\2f'^T^^^^F{X,)i® B{\,)F,{\,, A3) (C.8) 
pyAi — A2) 

— * 

which appears on the vector $3(Ai, A2, A3) and has to be canceled out by 

- K\.f'%^,F(M)i^ B(A3) 'S^^^^ (C.9) 
p(Ai-A2) /(A3-A2) 

and therefore we find that the function -F2(A2, A3) satisfies 

Hence, equations (C.4) and (C.IO) are able to fix the functions F2(A2, A3) and F3(A2, A3), and 
as a consequence follows the expressions (49) and (50) of section 5. We have also checked the 
consistency of these results by verifying that all other terms satisfy the condition of symmetry 
(51) of section 5 . Such calculation involves additional properties, since we need to turn twice 
the operators B{\i) . Here we list some extra identities which are extremely useful to prove 
the symmetry under Ai •(-^ A2, 

[i® B{y)].r^2{x) = B{y), [B{y) ® i].r23{x) = \^B{y) ® i (C.ll) 
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[B{y) ^].ruix) = B{y) ® ^+ rllix)^® B{y), [^0 B{y)].r23{x) = ^0 Biy) + rllix)Biy) ^ 

(C.12) 

Finally, another useful test is to look for certain unwanted terms which must be automati- 
cally canceled out. For instance, this is the case of the terms 

[/(A2)]^i?a(A)i?a(Ai)i?a(A3), [liX^)]"- Ea{X) E,i\,) B,{X,) (C.13) 

It is direct to verify that the first term in eq.(C.13) is canceled out with the help of F2{X2, A3) 
while the second term depends on the functional form of the function -F3(A2, A3) . 
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Figures 

Figure 1. The 36 nonvanishing Boltzmann weights of the rational spl{2\l) model. 
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